Raman laser spectroscopy of Wannier Stark states 
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Raman lasers are used as a spectroscopic probe of the state of atoms confined in a shallow 
ID vertical lattice. For long enough laser pulses, resolved transitions in the bottom band of the 
lattice between Wannier Stark states corresponding to neighboring wells are observed. Couplings 
between such states are measured as a function of the lattice laser intensity and compared to 
theoretical predictions, from which the lattice depth can be extracted. Limits to the linewidth of 
these transitions are investigated. Transitions to higher bands can also be induced, as well as between 
transverse states for tilted Raman beams. All these features allow for a precise characterization of 
the trapping potential and for an efficient control of the atomic external degrees of freedom. 



Cold atoms trapped in optical lattices have proven to 
be well suited for simulating solid-state systems, making 
possible the observation of Bloch oscillations [l|, , res- 
onant tunneling [H-IH or the Mott insulator regime Q. 
Besides, the precise knowledge and control of the atomic 
external degrees of freedom in these systems make them 
promising for applications such as metrology Q or 
the development of inertial sensors [1, [l^. In a recent 
article [ll| we showed that Raman pulses can be used 
to induce tunneling between neighboring sites of a ID 
vertical lattice. The present article aims at providing a 
more detailed description of the system. In particular we 
use Raman spectroscopy to probe the energy structure of 
atoms trapped in a ID vertical lattice. This study is mo- 
tivated by recent proposals to use such a system for short 
range forces measurements 1^ 15L or for the realization 



of a compact gravimeter [Tc , 16 , 17 1 . It is also of interest 
for any experiment using a shallow optical lattice, as it 
provides a comprehensive characterization of the system. 

We consider atoms trapped in a vertical standing wave, 
created by a laser far detuned from resonance. This re- 
sults in a periodic potential, which is superimposed to 
the gravitational potential in the vertical direction. The 
internal atomic structure is approximated by a two-level 
system with long lived states \g) and |e) with energy dif- 
ference hveg- The total Hamiltonian of this system is 
given by 



H — Hint + Hi + H„ 



(1) 



where Hint = hiyeg\e){e\ represents the internal energy, 
Hi = Uo{l — cos(2fc;f ))/2 is the periodic lattice potential 
with lattice depth Uq, lattice wavenumber ki and verti- 
cal spatial coordinate z, and Hg — niagz represents the 
gravitational potential, where rua is the mass of the atom 
and g is the gravity acceleration. 

As known from solid state physics, the eigenstates of 
the external part Hi + Hg, which is the sum of a peri- 
odic and a linear potential, are given by Wannier-Stark 



(WS) states [H, They form a so-called WS-ladder 
of states |Wb_m), where b is the discrete well index of 
the Bloch band, which structures the eigenstates in the 
periodic lattice, and the discrete quantum number m is 
the well index, characterizing the vertical position of the 
wave function {z\Wb.m) and which labels the well contain- 
ing its main peak in the limit of deep lattices (Uq ^ Er). 
The energy difference between adjacent lattice sites in 
the same band is simply the potential energy difference 
between two neighboring wells hvB — 'Tiag^i/'^, where 
A; = 271 /ki is the lattice wavelength, and vb is the Bloch 
frequency. Considering only the bottom band (6 — 0) 
and adding the internal Hamiltonian leads to a new WS- 
ladder like structure consisting of states \g, Wm) and 
\e,Wm') (see figure[T]). 
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FIG. 1: Wannier-Stark ladder of states and couplings between 
states by the probe laser 

In this structure, the application of a laser field reso- 
nant to A£'//i = Veg -I- Amj/B with Am = 0, ±1, ±2, . . . 
allows for coupling one state of this ladder to neighbor- 
ing WS states with opposite internal state, and thereby 
for the determination of the clock frequency Vf,g and 
the local gravity g in a spectroscopic measurement. In 
this, the inherent state labeling [20| gives us a tool for 
the precise measurement of the external state by inter- 
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nal state detection. Coupling of the ladder states be- 
comes apparent when adding a coupling Hamiltonian 
Hs — hflug^o cos{ujt — ksz)\e){g\ + H.c. to H, where 
^Uo=o is the Rabi frequency in absence of the lattice 
potential and kg — 2tt/ the coupling laser's wave num- 
ber [2l|. From this, the coupling strength for transi- 
tions between pairs of these states either in the same well 
(Am = m — m' = 0), or in neighboring wells (Am ^ 0) 
is calculated to be [2l| : 

nAm = nu„=o{Wm\e"'''\Wm'). (2) 

As we will see in more detail later in this paper, the lat- 
tice depth Uo plays an important role for driving these 
transitions. In too shallow lattices, the atomic localiza- 
tion is too weak and Landau-Zener (LZ) tunneling oc- 
curs, which limits the WS state lifetime. For too deep 
lattices, the WS states are localized in only one well, 
which strongly limits the intersite coupling strength and 
thus compromises spectroscopy measurements. Figure [2] 
illustrates the delocalization of the WS wavefunction at 
the depth of 1.6 Er we use in our experiment: it dis- 
plays the spatial density probability of WS wavefunction 
|^o(^)P = l(^l^o)Pi which extends over about 15 wells. 
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FIG. 2: Spatial density distribution of the \ Wo) wavefunction 
for a lattice depth of Uo = 1.6 iSr- 

Intersite transitions can be realized choosing As in the 
optical domain, with a single laser connecting two dif- 
ferent electronic states, such as the ground state and a 
metastable state of the optical clock transition studied in 
in . Alternatively, they can be realized with two pho- 
ton transitions between the two hyperfine ground states 
of alkali atoms. In this paper, we will focus on this latter 
case. For *^Rb, transitions between the ground and ex- 
cited hyperfine levels \g) — |5^S'i/2,F — l,m_F = O) and 
|e) — |5^S'i/2, = 2,mp — O) can be driven using coun- 
terpropagating vertical Raman beams providing a fre- 
quency difference of — 1^2 ~ i^i, that can be tuned 
around v^g = 6.834 GHz. This transition implies a mo- 
mentum transfer k = ki -\- k2 ~ 47r/As with = 780 nm 



replacing kg in equation[2] Here, Vi^2 and fci.2 are the re- 
spective frequencies and wavenumbers of the two Raman 
lasers. 

The intersite coupling on this transition is discussed 
for different values of A; in section ID The experimen- 
tal apparatus is then presented and the observed cou- 
pling strengths are compared to theoretical values. The 
achieved linewidth surpassing the Fourier limit by less 
than a factor two in the range up to 1.4 s of spectroscopic 
interrogation time and its limitations are presented and 
discussed in section [Hi Finally, we show the observed 
longitudinal and transverse structures observed in our 
composed trap in section iHll 



I. COUPLINGS 

The possibility to drive resolved intersite transitions 
strongly depends on the lattice wavelength and depth, 
as illustrated in figure [3] We calculated the coupling 
strengths as a function of the lattice depth (using a nu- 
merical calculation of the second term in equation 2, 
see also [2l|) for Am = m! — ni = ±1 transitions, for 
three different lattice wavelengths: close to resonance 
(A; — 800 nm), far blue-detuned (A/ = 532 nm) and 
far red-detuned (A/ — 1064 nm). The choice of these 
wavelengths is motivated by the availability of powerful 
enough lasers (> 20 W in the far detuned cases, up to 
several watts close to resonance) , allowing to reach a suf- 
ficient depth of a few recoil energies with a relatively large 
waist (w 1 mm). Such a waist is required for minimiz- 
ing parasitic forces due to the dipolar potential gradient 
along the longitudinal direction when using such a sys- 
tem for high precision measurements (gravimetry, short 
range forces, fine structure constant [23|...). 
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FIG. 3: Coupling strengths Am = ±1 versus lattice depth 
for three different lattices wavelengths A; = 532 nm, A; = 
1064 nm and A; = 800 nm. K^l&Q nm. 
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At sufficiently low depths, of a few Er, we find cou- 
pling strengths of the same order of magnitude for the 
far blue and red detuned cases whereas the intersite cou- 
pling remains small for the close to resonance case. We 
find relatively large variations and modulations for the 
far detuned cases, with respect to the close to resonance 
case. Calculations performed for Am > 2 transitions 
showed similar behaviors. Nevertheless, tilting by the 
same angle both Raman lasers with respect to verticality 
allows to preserve the direction of k while reducing its 
magnitude, thus changing the coupling strengths. As an 
example, we plot in figure 2] the couplings in the close 
to resonance case for k reduced by a factor of 2, which 
corresponds to an angle of 60° with respect to verticality. 
Although we find larger coupling strengths, which oscil- 
late versus lattice depths, couplings comparable to the 
far detuned cases for Am > 1 are only reached at twice 
lower lattice depth. 
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FIG. 4: Coupling strengths Am = 0, ±1, ±2, ±3 versus lattice 
depth for a lattice wavelength A; = 800 nm and tilted Raman 
beams giving an effective wavelength of As = 1600 nm. 

Another important parameter depending on the 
lattice wavelength is the tunneling rate out of the 
bottom band, which limits the lifetime of the atoms in 
shallow lattices. For an estimate of this rate, we use the 
Landau-Zener formula. Figure [S] displays the calculated 
rates as a function of the lattice depth (in units of recoil 
energies) for low depths of 2 and 3 Er . The LZ tunneling 
rate remains small in the far blue detuned case, even 
for lattices as shallow as 2Er, whereas it becomes 
comparable to 1 s^^ for a lattice depth Uq between 2 
and 3 Er in the close to resonance case. Operation at 
A; — 532 nm thus appears more appealing as one can 
combine large lifetimes and good couplings for large site 
offsets Am. 

We have compared these calculations with measure- 
ments corresponding to the blue detuned case. Our sys- 
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FIG. 5: Landau Zener tunneling rate as a function of the 
lattice wavelength for Uo ~ 2 and 3Er. 



tem consists in laser-cooled Kh atoms loaded in the 
first band of a vertical one-dimensional optical lattice, 
created by a single mode frequency doubled Nd:YV04 
laser (A/ = 532 nm, maximal power 12 W) with a waist 
of about 700 ^m. As this blue detuned standing wave 
does not provide transverse confinement, a red detuned 
(A = 1064 nm, beam waist 200 /im) Yb fiber laser is su- 
perimposed to the lattice (see figure ^ . The difference 
in the waists of the two lasers allows to reduce inhomo- 
geneities in the lattice depth due to the transverse ex- 
tension of the atomic sample. Before being transferred 
into this mixed dipole trap, about 10^ atoms are loaded 
in a 3D-Magneto-Optical trap (MOT) and cooled down 
to 2 fiK with a far detuned molasses. The dipole trap 
lasers are switched on either at the end of this cool- 
ing phase, or at the beginning of the MOT sequence. 
Switching off the molasses lasers leaves about 1 % of 
the atoms trapped in the mixed trap with a lifetime of 
about 1 s. These atoms initially distributed in all the Zee- 
man sublevels of \5'^Si/2,F = 2) are then depumped to 
|5^S'i/2,F — l) before being optically pumped (98% ef- 
ficiency) on the \5'^Si/2,F ^ l) ^ \5^P^/2,F ^ O) tran- 
sition to the 15^51/2, = l,mF = O) Zeeman sublevel, 
which is sensitive to stray magnetic fields only to sec- 
ond order. After being released from the optical trap, 
atoms fall for about 140 ms before reaching the detec- 
tion zone located at the bottom of the vacuum chamber. 
The detection scheme is based on a time of flight mea- 
surement similar to the ones used in atomic clocks and 
inertial sensors. It allows to measure by fluorescence the 
atomic populations in the two hyperfine states F = 1 and 
F = 2, denoted A^l and N2, respectively [l^], from which 
we derive the transition probability P = N2/{N2 + Nl). 
The Raman transitions are driven by two counterprop- 
agating, circularly polarized beams at 780 nm, detuned 
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FIG. 6: Experimental setup. The laser beams for optical trap- 
ping (lattice at 532 nm and transverse confinement at 1064 
nm) and Raman spectroscopy (780 nm) are superimposed us- 
ing dichroic optics. 



from the atomic transition by about —3 GHz, and aligned 
along the direction of the optical trap beams. These are 
coUimated with a 1/e^ radius of 1 cm, ensuring a good 
intensity homogeneity along the transverse size of the 
trap. 

In order to determine the coupling strengths, a Ra- 
man spectrum is first scanned by measuring the transi- 
tion probability as a function of the frequency difference 
between the Raman lasers v^. For such scans, the in- 
tensities in the Raman laser beams are 0.25 and 0.54 
mW/cm^, and the duration of the Raman pulse is 8 ms. 
This ratio between the Raman intensities is chosen to 
cancel the differential light shift they induce on the fre- 
quency of the hyperfine transition (25| . 

We observe multiple resonances, corresponding to 
transitions between the two hyperfine levels at Raman 
frequencies equal to the hyperfine splitting plus or minus 
an integer number Am of Bloch frequencies {vb ~ 569 Hz 
in our system). The difference in the peaks heights is 
due to the difference in the coupling strengths. We then 
fix the Raman frequency difference at the center of each 
peak, record a Rabi oscillation pattern by measuring the 
transition probability as a function of the pulse length. 
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FIG. 7: Raman spectrum showing the transition probability 
as a function of the Raman frequency from a lattice depth 
of 1.6 Er- The resonances separated by the Bloch frequency 
vb ~ 569 Hz are the signature of intersite transitions. 



from which we extract the Rabi frequency. We repeat 
this procedure for different lattice laser power values. 
The measured coupling strengths of the Am = 0,1,2,3 
transitions are plotted in figure [S] as a function of the 
lattice laser power. The results have been normalized 
and the relation between lattice laser power and actual 
lattice depth has been adjusted for a better match with 
the theoretical predictions. Data points corresponding to 
minima in the couplings have larger error bars as inho- 
mogeneities in the Rabi frequencies damp so heavily the 
Rabi oscillations that almost no oscillations are observed. 
We attribute this inhomogeneity to the transverse spread 
of the atoms, which experience different lattice depths 
due to the Gaussian profile of the lattice beam. 

The good agreement between the measurements and 
the theoretical predictions allow us to determine the 
depth with a resolution on the order of 0.1 Er from the 
direct comparison of the relative amplitudes of the peaks. 
Alternative techniques for the determination of the lat- 
tice depth are not as accurate here: diffraction in the 
thick grating limit ^] creates sidebands in the veloc- 
ity distribution which can hardly be resolved due to the 
width of the initial velocity distribution of the order of 
2.5vr, and parametric excitation gives rise to very wide 
resonances due to the complete anharmonicity of the lat- 
tice potential. 



II. LINEWIDTH 

We then investigated the question of the linewidth of 
such transitions. Various effects are expected to con- 
tribute to the broadening of the transitions and to ul- 
timately limit the minimally attainable linewidth. One 
of them is the differential light shift (DLS) induced by 
the trapping laser beams. This effect is dominated by 
the light field of the transverse dipole trap, as atoms are 
trapped at the intensity maxima of the 1064 nm beam 
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FIG. 8: Normalized Rabi frequencies measured as a function 
of the lattice depth, for Am = 1,2,3,4. The normalization 
factor f2(7o=o is an adjustable parameter. Solid lines are the 
result of numerical evaluation of equation [21 



and observe maximal DLS/i? at its center. By perform- 
ing microwave spectroscopy on the 1 5^51/2, F = \,mp — 
0) — >■ |5^S'i/2,-F — = 0) transition, we measure a 

shift in the center of the line of about 3 Hz/W, and a 
broadening of about 2 Hz/W, which gives a limit to the 
Hnewidth of 3 Hz at the 1.5 W we typically use. We have 
also measured the shift of the line induced by the lattice 
laser beam and found a much smaller effect of 0.4 Hz at 
full power. 



the evolution of the linewidth as a function of the du- 
ration of the microwave pulse, where the microwave 
power has been adjusted for optimal transfer at reso- 
nance (which corresponds to the case of a tt pulse), as 
well as in figure I10[ which shows as a dotted line the mi- 
crowave spectrum corresponding to a pulse duration of 
1.4 s. 
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FIG. 10: Transition probability as a function of the microwave 
detuning with (continuous line) or without (dotted line) the 
DLS compensation beam (see text), for a microwave pulse 
of T = 1.4 s. The compensated transition's detuning to the 
hyperfine frequency of 18.7 Hz is the Zeeman quadratic shift 
due to a bias field of 180 mG. The mean DLS imposed by the 
transverse trapping laser is —4.6 Hz. 




10 100 1000 

Pulse duration (ms) 



FIG. 9: Linewidth (FWHM) of transitions between the hy- 
perfine states driven by a microwave probe in the same lattice 
well or by counterpropagating Raman lasers with Am = 3, 
as a function of the pulse duration r. For each different pulse 
duration, the probe's Rabi frequency is adjusted so that 
Q3T = n. The solid line shows the theoretical Fourier-limited 
FWHM of a pulse of duration r. 

This broadening is illustrated in figure [HI which shows 



The DLS/fl induced by the transverse trapping laser 
can be compensated thanks to an additional laser beam 
with a blue detuning for the |e) state and a red detuning 
for the \g) state [27l|. For that purpose, a small fraction 
of one of the two Raman beams is used, with an addi- 
tional detuning of 80 MHz in order to prevent undesired 
Raman transitions. This beam is superimposed with the 
transverse trapping laser beam and its size, position and 
power are adjusted to reduce the broadening of the mi- 
crowave transition. For a transverse trapping laser power 
of 1.5 W, the differential light shift is compensated with 
a power of 12 nW. Figure [TU] displays as a continuous line 
the microwave spectrum for optimal compensation. 

The linewidths of the Raman transitions are displayed 
in figure |9l for uncompensated (resp. compensated) 
DLS/i^, as squares (resp. circles). The ratio between the 
intensities of the two Raman lasers is set to cancel (on 
average) the net differential light shift they induce. How- 
ever, due to differences in the spatial modes of the two 
Raman lasers and parasitic reflections, this compensation 
is not perfect, which leads to a broadening of the hyper- 
fine transition. As this broadening is proportional to the 
total Raman intensity, the linewidth is proportional to 
the Rabi frequency of the transition, only increased with 
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respect to the Fourier-limited microwave transition by a 
constant factor of about 1.5. Atom loss prevents us from 
driving longer transitions. Nevertheless, our system al- 
lows us to achieve a spectroscopic resolution of about 
1 Hz, which can be of interest for selecting atoms in a 
single site of the lattice, as demonstrated in [2I] with 
less resolved transitions. 

In order to study the short term sensitivity of our sys- 
tem, we performed a spectroscopic measurement of the 
Bloch frequency by measuring alternately the frequency 
of the Am = -|-3 and Am = —3 and calculating the differ- 
ence to cancel any shift of the hyperfine clock frequency. 
We obtained a statistical uncertainty on the Bloch fre- 
quency of 2 X 10^^ in relative value after 1 s of integra- 
tion, which is a factor of 3 better than the previously 
reported sensitivity in [11] using Ramsey spectroscopy. 
The best sensitivity reported for a trapped accelerome- 
ter was 1.4 x 10~^ in relative value after one hour mea- 
surement time jTo| . which corresponds to an equivalent 
relative short term sensitivity of 9 x 10~^ at 1 s. 



III. PROBING THE LATTICE'S 
LONGITUDINAL AND TRANSVERSE 
STRUCTURES 

Performing the Raman scan over a larger frequency 
range reveals additional features. Figure [TT] displays such 
a spectrum for two different lattice depths of 4 and 2.3 Er 
(measured by Raman spectroscopy as described in the 
last paragraph of section|I]) . Note the large blue sideband, 
which corresponds to transitions from the bottom band 
to the first excited band. The absence of red sideband 
indicates that the trapped atoms have been loaded in the 
bottom band and that upper bands are not populated. 

We have measured the lifetime of the atoms in the first 
excited band for a lattice depth oi 4Er. A Raman pulse 
of 2 ms detuned by 30 kHz transfers 60 % of the atoms 
initially in = 1 in the excited band in the F = 2 state. 
When increasing the delay between the Raman pulse and 
the turning off of the lattice laser, and measuring the 
number of the atoms that have remained trapped, we 
observe a decay in the number of atoms in = 2. Cor- 
responding data are displayed in figure I12[ from which 
we extract an exponential lifetime of 16 ms. This rela- 
tively short lifetime explains that for trapping times as 
large as several hundreds of ms, only the bottom band is 
populated. Loading the shallow lattice from the initial 
thermal distribution simply selects atoms loaded in the 
bottom band. 

Transitions between vibrational states along the trans- 
verse directions can also be induced by the Raman lasers, 
provided that their effective wavevector projection along 
the transverse direction is not null. Such transitions 
are exploited for instance for Raman sideband cooling 
[29I [soj . To do so, we simply tilt the Raman retrore- 
flecting mirror. Figure [T^] displays zooms of the Raman 
spectrum close to the Am = -f3 transition, for a Ra- 
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FIG. 11: Transition probability as a function of Raman fre- 
quency. The broad peak at vanishing relative Raman laser 
detuning {uRVeg = corresponds to unresolved intersite tran- 
sitions in the same lattice band. The large structure arising 
between 10 and 50 kHz is due to a coupling from the bottom 
to the first excited band. 
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FIG. 12: Population in the first excited band as a function of 
trapping time. The solid line is an exponential decay fit to 
the data from which we extract a lifetime of 16 ms. 



man pulse of 400 ms, for k perfectly vertical (continuous 
lines) and tilted by a few mrad (dotted line). We find 
red and blue sidebands about 25 Hz apart from the car- 
rier, which correspond to transitions An = ±1, where 
n is the index of the transverse vibrational level. We 
find equal amplitude for both sidebands, indicating that 
atoms are distributed among many such n states. The 
transverse temperature was independently measured by 
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time of flight fluorescence imaging to be 1 fiK. In addi- 
tion, the sidebands are significantly broadened with re- 
spect to the carrier. This broadening is attributed to the 
anharmonicity of the potential, the depth of the trans- 
verse dipole trap being only about 4 times the average 
transverse kinetic energy. 
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FIG. 13: Dotted line: Transition probability as a function 
of Raman frequency around the Am = -f 3 transition, when 
the Raman lasers wavevector is slightly misaligned from the 
lattice's one. The two sidepeaks correspond to intersite tran- 
sitions involving a change in the transverse vibrational state. 
Solid line: the Raman lasers wavevector is aligned with the 
lattice's one. 



IV. CONCLUSION 

We have investigated the possibility to drive intersite 
transitions in an optical lattice using Raman transitions. 
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